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Abstract 

We put independent model dynamical constraints on the net electric charge Q of some 
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j^j} astronomical and astrophysical objects by assuming that their exterior spacetimes are de- 

scribed by the Reissner-Nordstrom metric, which induces an additional potential C/rn ^ 
Q 2 r~ 2 . From the current bounds Avj on any anomalies in the secular perihelion rate m 
of Mercury and the Earth-mercury ranging Ap, we have \Qq\ < 1 — 0.4 x 10 18 C. Such 
constraints are ~ 60 — 200 times tighter than those recently inferred in literature. For 

■ the Earth, the perigee precession of the Moon, determined with the Lunar Laser Rang- 
ing (LLR) technique, and the intersatellite ranging Ap for the GRACE mission yield 
\Q®\ ^ 5 - 0.4 x 10 14 C The periastron rate of the double pulsar PSR J0737-3039A/B 

\ system allows to infer |Qns| 5 x 10 19 C. According to the perinigricon precession of the 

main sequence S2 star in Sgr A*, the electric charge carried by the compact object hosted 
in the Galactic Center may be as large as \Q,\ < 4 x 10 27 C. Our results extend to other 
hypothetical power-law interactions inducing extra-potentials C/ pcr t = 'f'r^ 2 as well. It turns 
out that the terrestrial GRACE mission yields the tightest constraint on the parameter iff, 

■ assumed as a universal constant, amounting to |$| < 5 x 10 9 m 4 s~ 2 . 

d ; PACS: 04.80.Cc, 95.10.Km, 96.30.-t, 91.10.Sp, 97.60.Jd, 97.60.Lf 



1 Introduction 

The Reissner-Nordstrom (RN hereafter) metric [TJ [2J E] is the unique spherically symmetric 
and asymptotically flat solution of the Einstein-Maxwell equations. It describes the exterior 
spacetime around an isolated spherical object of mass M and electric charge Q. In standard 
isotropic coordinates [3], commonly used for processing observations in several astronomical 
and astrophysical scenarios, the RN metric coefficients are 

500 = 1 + ^r, 

g 0i = 0, i = 1,2,3, (1) 

511 = 522 = 533 = - (1 - 7^) , 
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where c is the speed of light in vacuum, and U denotes the gravitational potential. For an 
isolated spherical body, it is 



U-U N + U RN - _ + _ 5 - J (2) 

where Q is measured in C, G = 6.67384 x 10~ n N m 2 kg -2 is the Newtonian constant of 
gravitation, and 

1 q N m 2 

k e = - = 8.98755 x 10 9 — =-, (3) 

47T£() t> 

in which e = 8.854187817 x 10~ 12 N _1 m~ 2 C 2 is the electric constant. Thus, the following 
constant factor can be defined 

^ = |g = 3.3363 xl0--^. (4) 

Concerning the physical relevance of C/rn i n astronomical and astrophysical contexts, it 
relies upon the existence of macroscopic bodies stably endowed with net electric charges [5]. 
Conversely, experimentally studying the orbital consequences of the Reissner-Nordstrom metric 
can allow, in principle, to measure, or, at least, constrain the value of the total net Q carried 
by macroscopic bodies of astronomical size in a dynamical independent model way. 

The question whether the Earth carries a net electric charge is quite old [BJE]. Nonetheless, 
it is rather poorly defined [8]. Indeed, it is often unclear if one refers to a quantity obtained 
by integrating over the solid and liquid Earth's surface-the "globe"-, or to the whole "planet" 
including the atmosphere and exosphere as well, with a deliberately set outer boundary [8]. 
Clearly, it is the second case that matters for our scopes, but it is difficult to have reliable 
estimates for it [8]. Concerning the Earth's surface, it should carry a net electrostatic negative 
charge of the order of El [10] |Q© rf | ~ 4 — 5.7 x 10 5 C because the electrostatic potential 
near the surface of the Earth is about 100 V [£l [10]. According to the author of Ref. [7J, 
tentative measurements of the net electric charge of the Earth as a planet would indicate that 
the electric field in the near-Earth space varies from ~ 0.1 mV m _1 to ~ 10 mV m , implying 
IQel < 50 C. 

As far as Sun-like, main sequence stars are concerned, actually, they should not be exactly 
neutral, as recognized in the early twenties of the 20th century \12 \ [T3 l 114 1 I16j. Indeed, 
the average velocity of free electrons in a plasma in thermodynamics equilibrium is higher than 
that of free protons, so that a larger number of electrons than protons on the stellar surface 
should attain the escape velocity, thus creating a net imbalance of positive charge [171 [THl [18] . 
As a result, if mp,m e ,q e are the mass of the proton, the mass of the electron and the charge 
of the electron, respectively, the expected total charge carried by a spherical, nonrotating star 
can be evaluated as [18] 

Q* = ^-^ M, = 77.03 C (5) 

Qe 



for M* = Mq. Such a figure is in agreement with the bound set by the authors of Ref. [19] 
according to whom the net positive electric charge carried by a Sun-like, nonrelativistic star 
cannot exceed ~ 100 C. See also Ref. [20]. On the other hand, in 1913 Hale [IT], in the 
framework of his attempts to detect the general magnetic field of the Sun, argued that the 
residual volume charge needed to account for it should have a negative sign. 

Net electric charges in extreme astrophysical objects like white dwarfs, neutron stars and 
black holes may have profound influences on several fundamental processes occurring in such 
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scenarios |21|,I22J . E.g., an overcritical electric field could occur in the region named dyadosphere 
[23] around Reissner-Nordstrom black holes, leading to the creation of e + — e~ pairs out of 
the vacuum which can cause a [23J Gamma Ray Burst (GRB) The authors of Ref. [SB] 

deal with nonrotating neutron stars in hydrostatic, chemical (beta), and diffusive equilibrium. 
Baryons have high masses and are, therefore, kept inside the star by its gravitational potential 
well, whose depth must be larger than the kinetic part of their Fermi energies. Anyway, such 
a potential is not deep enough to hold the notably lighter electrons, whose Fermi energy is 
relativistic (much larger than m c 2 ) and which would therefore tend to escape, yielding a net 
positive charge on the neutron star. However, at the typical densities occurring in neutron stars, 
already a small charge imbalance between protons and electrons causes an electrostatic field 
that can hold in the electrons and at the same time partially balance the gravitational force on 
the protons, preventing them from sinking into the center of the star, because of their relatively 
small kinetic energies. Thus, such models of neutron stars contain a spherically symmetric 
electrostatic potential evaluated as large as U c ~ 10 8 V corresponding to a charge of 

Qns ~ Hl-26 C (6) 

by assuming i?NS = 10 km. The charge separation can occur also in rotating neutron stars en- 
dowed with magnetic field and surrounded by plasma like pulsars [27] • Although the model by 
the authors of Ref. [27] presents shortcomings from the point of view of a realistic description 
[28], it is nonetheless useful to illustrate some basic properties. Due to the complex interaction 
between the rotation of the star, its magnetic field and the surrounding plasma in its magneto- 
sphere, a p c (r,9) charge density occurs, where 6 is the colatitude. Assuming complete charge 
separation, the corresponding average number density is [29j 



i 

P\ 2 



(n e ) = 1.75 x 10 16 m- a ( - ) 1-^1 , (7) 
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where P is the rotation period of the neutron star. For "normal" pulsars, i.e. for [29j P ~ 
0.5 s, P ~ 10" 15 , Eq. ([7D yields (r = R m = 10 km) 

Q NS ~ 1.7 x 10 10 C, (8) 

while for millisecond pulsars with [29] P ~ 3 ms, P ~ 10~ 20 it is 

Qns ~ 6.8 x 10 8 C. (9) 

For highly compact stars, whose radius is on the verge of forming an event horizon, the extremely 
high density and relativistic effects may, in principle, affect also the net charge of a compact star, 
so that it can take much more charge to be in equilibrium [301 ED E21 E3 EH EH] . The authors 
of Ref. [36] pursue the ideas of the aforementioned works by dealing with the effect of charge 
in cold compact stars, made of neutrons, protons and electrons They assume that the charge 
density is proportional to mass density, and that the net charge in the system is in the form of 
trapped charged particles carrying positive electric charge. By posing Qns ~ y/GineoM, the 
authors of Ref. [36] find Qns ~ 2 x 10 20 C. It turns out [36] that such highly electrified stars 
would not be stable, collapsing to a black hole: during the collapse, very little amount of charge 
would have the opportunity to escape, so that a charged black hole would form. Lensing by a 
charged neutron star, whose exterior spacetime is modeled by the Reissner-Nordstrom metric, 
was studied in Refs. [53| [33]. 
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The authors of Ref. [55] investigated charged frozen stars. They are quasiblack holes with 
pressure, i.e. objects whose boundary approaches their own gravitational radius as closely as one 
likes, without actually reaching it because of strong internal relativistic effects. Such spherically 
symmetric relativistic charged fluid distributions are bounded by a physical surface of radius 
Rq . Their internal region is filled with a fluid characterized by its mass and charge densities and 
by a nonzero pressure. The exterior spacetime is described by the Reissner-Nordstrom metric 
[55]. 

Electrically charged strange quark stars were studied in Ref. [56]. Such objects would be 
compact stars made of absolutely stable strange quark matter. They may carry ultrastrong 
electric fields on their surfaces, of the order of ~ 10 20 V m" 1 for ordinary strange matter. 
Under certain circumstances, the strength of the electric field may increase to values that 
exceed 10 21 V m -1 . 

For electric charges in black holes, see |37 |. I38 |, [39 ], 140) . The formation of Reissner-Nordstrom 
black holes was intensively studied in recent years, also because the astrophysical conditions 
leading to it may look rather problematic: see, e.g., Refs. gTJ SZl S31 S31 S3 [361 SSI S3 
SHI S3 EQ] EI] [52] ■ Anyway, the existence of such a kind of black holes is neither forbidden 
by theoretical nor observational arguments [57]- Also the possibility that electrically charged 
rotating black holes, described by the Kerr-Newman (KN) spacetime metric [581 159] . really 
exist in nature is somewhat controversial |60| . Gravitational lensing in Reissner-Nordstrom 
spacetimes was studied in Ref. [61] on the basis of the results of Refs. [62] [63] [64] [65] . For some 
features of Schwarzschild black hole lensing, see, e.g., Ref. \66\ 167], Methods to measure the 
electric charge of black holes through retrolensing with the currently ongoing [68] astrometric 
mission Spektr-R (known also as Radioastron) [69] were proposed in Ref. [57]; see also Refs. 
[70"| [71] , There is nowadays wide consensus that the compact object with M, = 4 x 10 6 M Q 
hosted by the Galactic Center in Sgr A* [72] is a supermassive black hole. It is orbited by 
several main sequence S stars: the closest one so far discovered, known as S2, revolves around 
Sgr A* at (r) = 1433 AU from it in 15.98 yr [73]. The author of Refs. [73] [75] pointed out 
that lensed images of the stars orbiting close to Sgr A* can provide insight into the form of its 
exterior spacetime metric and, in particular, on its electric charge Q, as well. Gravitational 
lensing of stars orbiting the supermassive black hole in the Galactic Center was studied in Ref. 
|76] as well. According to the author of Ref. [70] , for Sgr A*, a charged black hole could be 
distinguished from a Schwarzschild black hole with the space radio telescope Radioastron, at 
least if its charge is close to the extremal value 

Qf 1 = y/GAneoM. = 7.4 x 10 26 C (10) 

allowed to avoid a naked singularity. In general, black holes and naked singularities can be 
observationally differentiated through their gravitational lensing characteristics [75] [TjJ] [SO] . We 
recall that there are no logical arguments preventing existence of naked singularities; according 
to Ref. [H], this is an open question. In the case of a Kerr-Newman metric, Q^ xtr is modified 
by the black hole's spin J, as 

QT tr = \/G4ne Q (1 - xl)M. = (6.3 - 6.6) x 10 26 C (11) 

since for Sgr A* it is [77] 182] 

*- = ^~ a44 - a52 - (12) 

For a Reissner-Nordstrom black hole, its charge changes the size of the shadows up to 30% in 
the extreme charge case. Therefore, the charge of the black hole can be measured by observing 
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the shadow size, if the other black hole parameters are known with sufficient precision [70j . 
The emission and absorption properties of charged black holes were studied in Refs. |83|, [84"1 
1851 186]- Their characteristic pattern could be used as footprint of the black hole. For quantum 
mechanical uncertainties in measuring, among other things, the charge of a black hole, see Ref. 
|87j . The geodesies in the Reissner-Nordstrom spacetime were studied in Ref. [88] , 

The appearance of charged, traversable wormholes in the Reissner-Nordstrom metric is 
investigated by the authors of Ref. [89] . 

In Sec. [2] of this paper we will use well known and largely tested orbital motions (see Sec. 
12. ip of some natural (Sec. 12. 2j) and artificial (Sec. I2.3H bodies in the solar system to infer 
constraints on the net electric charges Q carried by the Sun and the Earth. In Sec. 12.41 we 
will also use the well known, and extensively studied, double pulsar binary system and the 
main sequence S2 star orbiting the compact object in Sgr A*. The resulting constraints are 
inferred in a phenomenological, model-independent way; the only assumption made is that 
the exterior spacetime of the sources of the gravitational field is described by the Reissner- 
Nordstrom metric. For the motion of electrically and magnetically charged particles in such 
a spacetime, see Refs. [9Q|, 191] . Finally, we stress that our results are not necessarily limited 
to the Reissner-Nordstrom spacetime, being valid also for other theoretical schemes yielding 
power-law interaction potentials U pcr t oc r~ 2 [92 . This point is treated in Sec. El 

2 Reissner-Nordstrom long-term orbital effects and compari- 
son with the observations 

2.1 Analytical calculation of the secular precession of the pericenter 

The long-period effects caused by C/rn on the motion of an electrically neutral particle of 
mass m can, e.g., be computed perturbatively by adopting the Lagrange equations for the 
variation of the osculating Keplerian orbital elements [93] : their validity has been confirmed 
in a variety of independent phenomena. Generally speaking, the Lagrange equations imply 
the use of a perturbing function 1Z which is the correction U per t to the standard Newtonian 
monopole i/jj oc r _1 . In the case U pcv t = £7rni the average over one orbital revolution of the 
perturbing function 1Z is straightforwardly obtained by using the true anomaly / as fast variable 
of integration. It is 

(K) = hGQ2 r , (13) 

2 c 2 a 2(!_ e 2)2 

where a is the semimajor axis and e is the eccentricity of the test particle's orbit. From Eq. 
(|13p and the Lagrange equation for variation of the longitude of pericenter w 



dw \ 1 



T-eM 



dt I n^a 2 

it turns out that w experiences a secular precession given by 

'dw\ k e GQ 2 k e G^Q 



am + ^m T om h (14) 



dt j 2c 2 n h a A (1 - e 2 ) 2c 2 m\ o§ (1 - e 2 ) 



(15) 



In Eq. (I14p -Eq. (|15|) . = \f GM / a 3 is the Keplerian mean motion; I, entering Eq. (I14|) . is 
the inclination of the orbital plane to the reference {x, y} plane. The longitude of pericenter 
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w = £1 + uj is a "dogleg" angle since it is the sum of the longitude of the ascending node 
which is an angle in the reference {x, y} plane from a reference x direction to the intersection 
of the orbital plane with the plane itself (the line of the nodes), and of the argument 

of pericenter uj, which is an angle counted in the orbital plane from the line of the nodes to 
the point of closest approach, usually dubbed pericenter. The precession of Eq. (|15p . which 
is an exact result in the sense that the assumption e ~ was not made, agrees with the one 
obtained by Adkins and McDonnel in Ref. [95] with a more cumbersome calculation. To 
facilitate a comparison between such two results, we note that, in general, the authors of Ref. 
[95] work out the pericenter advance per orbit A6 p : it corresponds to our Aw = (w) Pb, where 
Pb = 27r/nb is the orbital period. Moreover, in the potential energy V(r) = a_y_ ( _i)r - ^ +1 ^ 
of Ref. [95] it must be posed j = 1 and a_2 — > 2~ l c~ 2 k e GmQ 2 , while in A9 p (—(j + 1)) of 
Eq. (38) in Ref. [95] it must be set L — > a (I — e 2 ), Xi( e ) = 2. With such replacements, 
it can be shown that the advance per orbit of Eq. (38) in Ref. [95] corresponds just to our 
precession in Eq. (|15|) . The precession of Eq. (|15|) yields a pericenter advance per orbit Aw 
which, by posing L -> o(l - e 2 ), agrees also with the shift calculated by the authors of Ref. 
|96] in their Eq. (16) from the geodesic equations of motion, and in their Eq. (37) with the 
Laplace- Runge-Lenz vector }97[ 198] , The rate of Eq. (|15p also agrees with that obtainable 
from the fractional shift computed by Lemmon and Mondragon [99] as (2ir)~ 1 A(p = —€q, 
where e Q = 2e M r 2 Q r~ 2 , r M = 2GMc~ 2 , r 2 Q = k^Gc^, e M = GM C - 2 a" 1 (l - e 2 )" 1 . For 
earlier derivations of the pericenter precession in the Reissner-Nordstrom metric, see Refs. 
[TOOl [T0T1 [T02l[T03] . 

The analytical result of Eq. ([15 j) is confirmed by a numerical integration of the equations of 
motion for Mercury with, say, Qq = 10 19 C displayed in Figure [TJ both yield —24 milliarcsec- 
onds per century (mas cty~ x in the following). The choice of the numerical value adopted for 
the electric charge of the Sun is purely arbitrary, being motivated only by the need of dealing 
with easily manageable figures on the vertical axis. 

The authors of Ref. |104| computed the correction Pq to the Keplerian orbital period Pb 
due to the Reissner-Nordstrom metric in the framework of the studies about general relativistic 
effects on bound orbits of solar sails [105} 1106] . 

2.2 Constraints on Q Q from solar system planetary orbital motions 

The corrections Aw to the standard Newtonian-Einsteinian secular precessions of the longi- 
tudes of the perihelia are routinely used by independent teams of astronomers [107} 1108] as 
a quantitative measure of the maximum size of any putative anomalous effect allowed by the 
currently adopted mathematical models of the standard solar system dynamics fitted to the 
available planetary observations. Thus, Aw can be used to put constraints on the parameters 
like Q entering the models one is interested in. From Eq. (|15h it turns out that the tightest 
constraints come from Mercury, which is the innermost planet with a = 0.38 AU. 

Fienga et al. [108] . who used also a few data from the three flybys of MESSENGER in 
2008-2009, released an uncertainty of 0.6 mas cty -1 for the perihelion precession of Mercury, 
so that it is 

Q Q < 1.5 x 10 18 C. (16) 

The uncertainty in the pre-MESSENGER Mercury's perihelion extra-precession by Pitjeva [107] 
is about one order of magnitude larger (5 mas cty" 1 ). 

Tighter constraints on Qq come from the interplanetary Earth-Mercury ranging. Indeed, 
according to Table 1 of Ref. [108] , the standard deviation a^ p of the Mercury range residuals 
Ap, obtained with the INPOPlOa ephemerides and including also three Mercury MESSENGER 
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RN shift of Mercury's longitude of perihelion (Qs U n = 10 Q 




20 40 60 80 100 

t(yr) 



Figure 1: Centennial shift Aw, in milliarcseconds (mas), of the longitude of perihelion of 
Mercury due to the Reissner-Nordstrom potential of the Sun for, say, Qq = 10 19 C, chosen 
just for illustrative purposes, computed from the difference between two numerical integrations 
of the equations of motion in cartesian coordinates performed with and without C/rn- Both 
the integrations share the same initial conditions, retrieved from the NASA WEB interface 
HORIZONS at |http:/ /ssd.jpl.nasa.gov/horizons.cgij for the epoch J2000.0. It agrees with the 
analytical result of Eq. (fT5l) . 



flybys in 2008-2009, is as large as 1.9 m. A numerical integration of the equations of motion of 
Mercury and the Earth with 

Q Q < 4.2 x 10 17 C (17) 

yields a Reissner-Nordstrom range signal with the same standard deviation; Figure [2] displays 
it. 

Our model- independent, dynamical bounds for the electric charge of the Sun are ~ 60 — 200 
times more stringent than Q@ < 8.91 x 10 19 C inferred by the authors of Ref. [96]. 

2.3 Constraints on from the Moon and the GRACE spacecraft orbiting 
the Earth 

Remaining within the solar system, let us, now, constrain the electric charge of the Earth with 
natural and artificial bodies. 

The orbit of the Moon is accurately reconstructed with the Lunar Laser Ranging (LLR) 
technique |109j since 1969; Figure B-l of Ref. [110] shows that the residuals of the Earth-Moon 
range are at a cm-level since about 1990. The secular precession of the lunar perigee is known 
with an accuracy of about 0.1 mas yr _1 11121 IH3| II 14] ; thus, Eq. (|15p yields 

Q e < 5.1 x 10 14 C. (18) 

The Gravity Recovery and Climate Experiment (GRACE) mission [115] . jointly launched in 
March 2002 by NASA and the German Space Agency (DLR) to map the terrestrial gravitational 
field with an unprecedented accuracy, consists of a tandem of two spacecrafts moving along 
low-altitude, nearly polar orbits continuously linked by a Satellite to Satellite Tracking (SST) 
microwave K-band ranging (KBR) system accurate to better than 10 /jm (biased range p) [116| 
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Earth-Mercury RN range perturbation (Q Sun = 4.2 10 17 C) 




t(yr) 



Figure 2: 2008-2009 Earth-Mercury range shift Ap, in m, due to the unmodeled Reissner- 
Nordstrom potential of the Sun for Qq = 4.2 x 10 17 C computed from the difference between two 
numerical integrations of the equations of motion in cartesian coordinates performed with and 
without C/rn- Both the integrations share the same initial conditions retrieved from the NASA 
WEB interface HORIZONS at http://ssd.jpl.nasa.gov/horizons.cgi, for the epoch J2000.0. It 
is (Ap) = 0.04 m, <ta p = 1.9 m. 



and 1 fim s _1 (range-rate p) [1161 I117j . Studies for a follow-on of GRACE |118] show that 
the use of a interferometric laser ranging system may push the accuracy in the range-rate to 
a ~ 0.6 nm s _1 level. A numerical integration of the equations of motion for GRACE A/B, 
including also the mismodelled signal of the first nine zonal harmonics of geopotential [93] 
according to the global Earth gravity model GOCO01S |119j . shows that the SST range is more 
effective than the SST range-rate in constraining the Earth's electric charge for which it holds 

Q e < 4 x 10 13 C, (19) 

which is about one order of magnitude better than the lunar constraint. Figure [3] depicts the 
numerically integrated GRACE SST range signal due to C/rn and the aforementioned mismod- 
eled zonals. It may be interesting to notice that the GRACE-based dynamically inferred upper 
bound of Eq. (fT9|) almost corresponds to the net electric charge (Q® = 2.9 x 10 13 C, from 
U e = 4.1 x 10 16 V) that the Earth's surface should have carried to account for the terrestrial 
magnetic field in the early theory propounded in 1879 by Perry and Ayrton |120j . It was soon 
rejected by Rowland [121] because of the manifest lacking of the enormous effects implied by 
it. 



2.4 Constraints from the double pulsar and Sgr A* 

The periastron of the double pulsar PSR J0737-3039A/B |122| I123j can be used to constrain 
QnSi at least for this particular system. The present-day accuracy in measuring the secular 
precession of the periastron is 6.8 x 10~ 4 degree per year (deg yr _1 in the following) |124j : thus, 
a straightforward application of it to Eq. (I15D would give 

Qns < 6.9 x 10 18 C. (20) 
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GRACE RN SST range perturbation (Q Earth = 4 10 13 C) 
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Figure 3: Daily GRACE SST range shift Ap, in pm, due to the unmodeled Reissner- 
Nordstrom potential of the Earth (Q® = 4 x 10 13 C) and the first nine mismod- 
elled zonals of geopotential according to GOCO01S |119| computed from the difference 
between two numerical integrations of the equations of motion in cartesian coordinates 
performed with and without C/rn- Both the integrations share the same initial con- 
ditions cont ained in the files GNV1B_2003-Q9-14_A_00 and GNV1B_2003-Q9-14_B_00 re- 
trieved from |ftp: / / cddis.gsfc.nasa.gov/pub / sir / predicts / current / graceA_irvs_ 081202_0.gfz and 
ftp://cddis.gsfc.nasa.gov/pub/slr/predicts/current/graceB_irvs_081201_l.gfz. The epoch is 13 
September 2003. (See ftp:/ /podaac.jpl. nasa.gov/pub/grace/doc/Handbook_lB_vl. 3.pdf for the 
explanation of the GPS Navigation Data Format Record (GNV1B) format). It is (Ap) = 
—2 fxm, o"Ap = 2 /im. 
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Actually, the larger uncertainty in the theoretical expression of the general relativistic 1PN 
periastron precession must be taken into account as well. It is as large as 0.03 deg yr _1 |126j . 
so that it yields 

Qns < 4.7 x 10 19 C. (21) 

Notice that the constraint of Eq. (|2ip is smaller by almost one order of magnitude than 
VG47re M t ot = 4.4 x 10 20 C for the pulsar PSR J0737-3039A/B system. For other electromag- 
netic/gravitational effects occurring in binary systems hosting compact objects, see, e.g., Ref. 

The perinigricon of the S2 star, orbiting in 15.98 yr the compact object hosted in Sgr A* 
with M = 4x 10 6 M Q [72] at (r) = 1433 AU from it [73] can be used to put dynamical constraints 
on the electric charge of the attractive center. Indeed, an accuracy of [73] 0.05 deg yr" 1 can be 
inferred for the perinigricon precession of S2: thus, 

Q. < 3.6 x 10 27 C. (22) 

By propagating the uncertainties in the estimated parameters of the Sgr A*-S2 system [73] 
entering Eq. (|15p . it is possible to infer 

a Qm = 4.6 x 10 26 C. (23) 

In view of Eq. (|10j) . this implies that the present-day knowledge of the S2 orbital motion still 
leaves room for the possibility that the massive object in the Galactic Center is not a Reissner- 
Nordstrom black hole since its electric charge might be larger than the extremal value to avoid 
a naked singularity. Discovering and monitoring other stars, closer to Sgr A*, will be crucial to 
shed further light on such an issue. For other potentially competing orbital effects which may 
take place in dense stellar cluster at the Galactic Center, see. Ref. [127] . 



3 Constraints on a generic r 2 extra-potential 

The bounds previously obtained in terms of electric charges in the Reissner-Nordstrom space- 
time can also be used to constrain the parameter ^ of a generic r -2 extra-potential [92] 

£W = J, (24) 

provided that the replacement 

ieQ 2 -> * (25) 

is made: = L 4 T -2 . Here we will adopt a purely phenomenological approach by considering 
\P just as an universal constant independent of the body which acts as source of the (modified) 
potential; for a discussion of power-law potentials, see Ref. [92]. In a quantum field theory 
framework, such kind of interactions arise from higher-order exchange processes with simulta- 
neous exchange of multiple massless bosons. In particular, potentials oc r~ 2 are generated by 
the simultaneous exchange of two massless scalar bosons [128J. 

It turns out that GRACE yields the tightest constrain on considered as a universal 
parameter, amounting to 

|^| < 5.3 x 10 9 m 4 s" 2 . (26) 

Among the solar planets, the most stringent bound comes from the Earth-Mercury ranging 
with 

|^| < 6 x 10 17 m 4 s~ 2 . (27) 
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The double pulsar and the Sgr A*-S2 system yields |^| < 7.4 x 10 21 m 4 s~ 2 and |^| < 4.3 x 
10 m 4 s -2 , respectively. 

4 Summary and conclusions 

We analytically worked out the long-term precession of the pericenter of a test particle orbiting 
a central body of mass M and electric charge Q in the framework of the Reissner-Nordstrom 
metric. We obtained a non- vanishing, secular precession whose expression is valid for any value 
of the eccentricity e of the particle's orbit. 

Then, we compared our analytical result to latest observations in several astronomical and 
astrophysical scenarios to infer upper bounds on the electric charge of various objects acting 
as sources of the Reissner-Nordstrom field. As far as our solar system is concerned, we get 
Q Q < 4 x 10 17 C from Mercury's motion. The GRACE mission around the Earth yields 
Q® ^ 4 x 10 13 C. Planned follow-on of GRACE and spacecraft-based missions to Mercury 
like MESSENGER (ongoing) and BepiColombo (to be launched in next years) will allow to 
strengthen such constraints. The periastron advance of the double pulsar PSR J0737-3039A/B, 
corrected for the uncertainty in the standard 1PN precession, yields Qpsr < 5 x 10 19 C. The 
orbital motion of the main sequence S2 star around the supermassive black hole in Sgr A*, 
monitored during a time span larger than a full orbital period, allows to infer Q, < 4 x 10 27 
C. This bound could become tighter in future if stars orbiting Sgr A* faster than S2 will be 
discovered and monitored. 

Our results are not necessarily limited to the Reissner-Norstrom spacetime due to a charged 
body, being, instead, valid for any perturbing potential U per t oc r -2 through a universal param- 
eter ^ having dimensions of = L 4 T~ 2 . In this respect, the tightest constraints come from 
the Earth-GRACE system with |^| < 5 x 10 9 m 4 s" 2 . 
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